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ABSTRACT 


A  highly  linear  passive  floating  element  wall  shear  trans- 
ducer has  been  developed  which  uses  a  sensitive  frequency  modula- 
tion technique  to  measure  a  capacitance  change  proportional  to 
the  displacement  of  the  floating  element  for  small  deflections. 
The  floating  element  is  4  mm  square  surrounded  by  a  gap  nominally 
0.18  mm  and  protrudes  into  the  boundary  layer  9  ym.  Scaled  on 
inner  parameters,  each  side  of  the  sensing  head  measures  about 
290,  the  mean  gap  width  is  13,  and  the  element  protrusion  is  0.7. 
The  transducer  was  tested  in  a  turbulent  boundary  layer  on  a  flat 
plate  at  momentum  thickness  Reynolds  numbers  of  3140,  3500,  and 
3760  and  the  results  compared  to  wall  shear  stress  estimates  ob- 
tained from  the  law  of  the  wall.  The  transducer  was  found  to 
predict  the  steady  mean  wall  shear  stress  within  the  limits  of 
experimental  uncertainty  in  all  cases;  transducer  measurements 
of  mean  wall  shear  stress  for  momentum  thickness  Reynolds  numbers 
of  3140  and  3500  were  within  2  percent  of  the  actual  value. 
Transducer  measurement  of  the  unsteady  mean  wall  shear  stress  was 
severely  limited  by  a  mechanical  resonance  in  the  transducer  at 
63  Hz. 


The  operating  principles  and  general  characteristics  of  the 
transducer  are  discussed,  its  performance  is  evaluated,  and  re- 
commendations are  made  concerning  future  work. 
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NOMENCLATURE 


A     surface  area  of  the  floating  element; 


C     capacitance  of  the  cantilever-stationary  plate  pair 


C  friction  coefficient 


d     mean  spacing  between  stationary  and  cantilever  plates 
with  the  cantilever  undeflected 


e     uncertainty  or  error  in  the  quantity  C 


EI     the  flexural  rigidity  of  the  cantilever  beam 

F     generally,  a  force;   in  Appendix  C,   a   functional  ex- 
pression for  the  quantity  being  measured 


*      frequency  in  Hz;   in  Appendix  A,   a  force  distributed 
the  length  of  the  beam  element  (f(x,t)) 


fQ     center  frequency  of  astable  multivibrator 


f,     fundamental  mechanical  resonant  frequency  of   the   beam 
element 


G     gain  in  the  frequency-to-voltage  converter  and  instru- 
mentation amplifier,  V/Hz 


J     mass  moment  of  inertia  of  the  floating  element 

k      surface  roughness 

k      static  spring  constant  of  cantilever  beam  to  transverse 


s 


loads,  3EI/L 


3 


length  of  the  cantilever  plate  and  beam  element 


M     mass  of  the  floating  element  (assumed   to  be  a  point 
mass) 


m 


U 

00 


mass  of  the  beam 


q      the  total  head,  —  p  u2  '>      also  the  total  error   (Appen- 
dix  C)  2   °° 


R      electrical  resistance 


Re     Reynolds  number  based  on  x,   the  distance   from  the 
plate's  leading  edge 


Re     momentum  thickness  Reynolds  number 


sensitivity  of  the  transducer  in  either  V/Pa  or  V/ym; 
also  the  spacing  between  capacitor  plates  (Appendix  A) 


t  thickness  of  the  cantilever  plate  ;  also  time 

u      longitudinal  fluctuating  fluid  velocity 
free  stream  velocity 


TT  -      W 

UT     friction  veloci 


:lty,  l[-f- 


Vq      transducer  offset  voltage 
V(t)    transducer  output  voltage 


w      width  of  the  cantilever  plate;   also  (Appendix  A)   the 
weight  per  unit  length  of  the  beam 


W 


weight  of  the  floating  element 
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x  referring  to  the  flow,  the  streamwise  distance  from  the 
plate's  leading  edge;  referring  to  the  beam  element, 
the  distance  along  the  beam  from  the  base  of  the  beam 


referring  to  the  flow,  the  distance  normal  to  the  wall 
measured  from  the  wall;  referring  to  the  beam  element, 
the  deflection  of  the  beam  element  from  its  undeflected 
position 


z      the  deflection  of  the  beam  with  respect  to  a  stationary 
observer 

S      boundary  layer  thickness;   also  (Appendix  A)  the  varia- 
tional operator  used  in  the  variational  calculus 


6* 


X* 


displacement  thickness 


transverse  deflection  of  the  free  end  of  the 
cantilever;  also  the  percent  most  probable  error  (Ap- 
pendix C) 


v 
the  viscous  length  scale, 


U 

T 


n      beam  deflection  normalized  on  mean  capacitor  spacing  d 


o 


normalized  coordinate  along  beam,  x/L 
)      momentum  thickness 

1      transducer  inclination  angle  from  the  vertical 
>      floating  element  angle  of  rotation 


T„     wall  shear  stress;  mean  wall  shear  stress  refers  to  a 
w 

spatial  average  as  opposed  to  steady  wall  shear  stress 
which  implies  a  time  average 


v 


kinematic  viscosity 
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Superscripts 


identifies  the  quantity  as  being  nondimensionalized  in 
inner  coordinates;  velocities  are  normalized  on  U  and 
lengths  are  normalized  on  X*  T 


Subscripts 
00      refers  to  free  stream  properties 
w      refers  to  quantities  measured  at  the  wall 
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1 .   INTRODUCTION 


1.1  BACKGROUND 

For  many  centuries,  men  have  studied  the  force  of  drag  which 
resists  the  motion  of  a  body  through  a  fluid.  Todd  (1967)  re- 
ports that  as  early  as  the  late  fifteenth  century,  Leonardo  da 
Vinci  conducted  tests  on  three  ship  models  with  different  hull 
forms,  presumably  to  identify  the  most  efficient  design. 

The  dawn  of  the  "age  of  steam"  in  the  late  nineteenth 
century  motivated  researchers  to  better  understand  fluid  drag  in 
order  to  both  predict  and  minimize  its  effect.  As  sailing  ships 
were  supplanted  by  larger  and  faster  self-propelled  ships,  fric- 
tion drag  was  recognized  as  a  sizeable  contributor  to  the  overall 
fluid  drag.  William  Froude  (1872)  conducted  a  series  of  friction 
drag  measurements  in  his  towing  tank  at  Torquay,  England,  to  ana- 
lyze the  effect  of  friction  drag  on  ships.  Froude  built  planes 
from  wooden  planks  and  towed  them  edge-first  through  the  water, 
measuring  their  drag  with  a  dynamometer.  In  so  doing,  Froude  ac- 
tually determined  the  mean  wall  shear  stress  acting  over  the  sur- 
face of  each  plane. 

Since  1904,  when  Prandtl  presented  his  boundary  layer  theory 
to   the  Mathematical  Congress  in  Heidelberg,  understanding  of  the 
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wall  shear  stress  has  grown  considerably  in  sophistication. 
While  Froude  intuitively  grasped  that  the  distribution  of  wall 
shear  stress  over  a  body  in  a  fluid  flow  was  not  uniform,  theory 
now  predicts  that  the  wall  shear  stress  is  a  function  of  both 
space  and  time.  This  dynamic,  nonhomogeneous  nature  of  the  wall 
shear  stress  is  particularly  evident  in  the  transition  and  turbu- 
lent regions  of  the  boundary  layer. 

The  wall  shear  stress  is  also  an  important  boundary  layer 
parameter  in  its  own  right  since  it  is  a  characteristic  of  the 
inner  region  of  the  boundary  layer  (near  the  wall)  where  viscous 
forces  predominate.  Many  researchers  have  found  that  boundary 
layer  processes  scale  on  the  wall  shear  stress,  particularly 
those  related  to  turbulence  in  the  transition  and  fully  turbulent 
regions  of  the  boundary  layer.  One  example  is  the  familiar  "law 
of  the  wall"  in  which  the  turbulent  mean  velocity  profile  is  non- 
dimensionalized  on  mean  wall  shear  stress  and  kinematic  viscosity 
[see  Coles  (1968)  and  White  (1974)].  Burton  (1974)  found  that 
turbulent  wall  pressure  spectra  and  corresponding  fluctuating 
longitudinal  velocity  spectra  scaled  on  the  wall  shear  stress. 
Gedney  (1979)  was  unable  to  collapse  his  transition  wall  pressure 
spectra  because  of  inaccurate  determination  of  the  wall  shear 
stress  upon  which  they  were  nondimensionalized.  In  short,  ongo- 
ing fundamental  boundary  layer  research  indicates  that  more  pre- 
cise measures  of  the  wall  shear  stress  than  currently  available 
are  increasingly  necessary. 


1.2  WALL  SHEAR  STRESS 

Wall  shear  stress,  like  any  stress,  is  a  derived  quantity 
and  thus  cannot  be  measured  directly.  The  wall  shear  stress  is  a 
dynamic  quantity  defined  at  a  point  by  the  following  limit: 

F 
_  lim   _A 

Tw  "  A->0    A  (1.1) 

where  F  is  the  friction  force  acting  over  area  A.  Unless  the 
measurement  of  the  wall  shear  force  is  made  at  a  single  point, 
the  best  we  can  expect  to  measure  is  a  mean  wall  shear  stress 
averaged  over  some  measurement  area.  Measurement  techniques  used 
to  determine  wall  shear  stress  are,  in  reality,  merely  techniques 
which  measure  quantities  related  to  the  wall  shear  stress  from 
which  the  wall  shear  stress  is  inferred.  As  the  measurement  area 
approaches  a  point,  the  measured  wall  shear  stress  approaches  the 
true  wall  shear  stress  according  to  equation  1.1.  Small  measure- 
ment areas  thus  strengthen  the  inference  that  the  measured  shear 
stress  is  the  true  shear  stress,  underscoring  the  need  for  good 
spatial  resolution  if  the  true  shear  stress  at  a  point  is  the 
subject  of  investigation. 
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1.3  MEASUREMENT  TECHNIQUES 

Table  I,  adapted  from  Winter  (1977),  provides  a  brief  list- 
ing of  the  mean  wall  shear  stress  measurement  techniques  current- 
ly in  use.  No  single  technique  is  universally  acceptable  because 
each  has  restrictions  on  its  application  to  different  flow  condi- 
tions and  possesses  limited  versatility.  Further,  no  technique 
can  make  an  adequate  measure  of  the  true  wall  shear  stress,  owing 
to  excessive  spatial  averaging  or  inadequate  bandwidth.  The  hot 
film  probe,  using  the  heat  transfer  analogy,  is  the  only  device 
with  sufficient  bandwidth  to  measure  some  of  the  unsteady  char- 
acter of  the  mean  wall  shear  stress.  However,  its  spatial  reso- 
lution and  sensitivity  are  both  poor. 


TABLE  I 
Mean  Wall  Shear  Stress  Measurement  Techniques 


Analogies 
Heat  Transfer 
Mass  Transfer 


Liquid  Films 


Flow  About  Obstacles 
Preston  Tubes 
Stanton  Tubes 
Razor  Blades 
Steps  &  Fences 

Static  Pressure  Holes 


Floating  Element 


Momentum  Balance 


Velocity  Profiles  (Wall  Similarity) 


In  an  effort  to  improve  the  resolution  of  dynamic  wall  shear 
stress  measurements,  a  new  type  of  floating  element  transducer 
has  been  developed  and  tested.  This  device  uses  a 
capacitance-measuring  technique  of  exceptional  sensitivity  [first 
developed  by  Moeller  (1983)]  to  measure  the  deflection  of  a  small 
floating  element  in  the  plane  of  the  wall.  The  sensitivity  of 
the  transducer  permits  its  operation  with  such  small  motions  of 
the  floating  element  that  the  device  exhibits  excellent  lineari- 
ty. Thus,  the  powerful  techniques  of  linear  systems  analysis  may 
be  used  to  uniquely  and  simply  relate  statistics  of  the  transduc- 
er output  to  corresponding  statistics  of  the  mean  wall  shear 
stress. 


2.   TRANSDUCER  DESIGN 


2.1  GENERAL  DESCRIPTION 

The  transducer  is  a  passive  floating  element  instrument 
which  generates  a  voltage  proportional  to  the  displacement  of  the 
element  and,  indirectly,  to  the  shear  force  acting  on  the  ele- 
ment. This  is  unlike  most  floating  element  devices,  which  gener- 
ally employ  an  active  force-balance  principle  so  that  the  element 
remains  nearly  stationary.  The  work,  of  Kistler  [Paros,  (1970)], 
Allen  (1980),  Frei  and  Thomann  (1980),  Escudier,  et.  al.  (1982), 
and  Nguyen,  Dickenson,  et.  al.  (1984)  are  all  examples  of  the 
active  technique.  In  the  current  design,  sensitivity  has  been 
improved  to  the  point  that  the  floating  element  deflects  less 
than  5  u  m.   For  a  streamwise  deflection  of  5  y  m,   the   floating 

element  protrudes  only  0.63  y  m  into  the  boundary  layer,  which  is 

-2 
a  y+  of  5x10    .   The  protrusion  due  to  floating  element  deflec- 
tion is  thus  well  within  the  viscous  sublayer  and  can  be  neglect- 
ed. 

Figure  2.1  presents  a  schematic  of  the  transducer.  The  beam 
element  is  the  moving  element  of  the  transducer,  composed  of  the 
floating  element  and  a  thin  cantilever  plate.  The  floating  ele- 
ment, mounted  in  the  plane  of  the  wall  atop  the  cantilever  plate, 
is  sensitive  to  motion  in  only  one  direction.    Shear   forces   on 
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the  surface  of  the  element  in  the  sensitive  direction  deflect  the 
cantilever  plate,  producing  a  change  in  capacitance  with  respect 
to  a  second  plate  parallel  to  and  mounted  a  short  distance  from 
the  cantilever  plate.  This  change  in  capacitance,  proportional 
to  floating  element  deflection  for  small  displacements,  changes 
the  frequency  of  a  square-wave  oscillator.  The  transducer  thus 
produces  a  frequency-modulated  signal  proportional  to  the  float- 
ing element  displacement  which  is  subsequently  demodulated  and 
filtered  by  signal  conditioning  electronics. 


FLOW 


WALL— 7  _ 


FLOATING 
ELEMENT -J 

CANTILEVER 
PLATE-* 


1 


GAP 


STATIONARY 
PLATE 


^m^ 


FIGURE  2.1   Transducer  Schematic  -  General 
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2.1.1  Transducer  Mechanical 

The  transducer  is  assembled  within  a  body  of 
phenolic-impregnated  linen  as  shown  in  Figure  2.2.  A  shallow 
slot,  milled  in  the  body  halves  before  assembly,  contains  a  sta- 
tionary capacitor  plate  of  copper-clad  laminate  epoxied  in  place. 


FLOATING   ELEMENT 
(4  mm  x  4  mm) 


CANTILEVER 


ACCELEROMETER 
LOCATION 


-STATIONARY 
PLATE 
(15  mm  x  30  mm) 


LINEN    PHENOLIC 
BODY  (2   HALVES) 


FIGURE  2.2 


Transducer  Schematic 


The  movable  capacitor  plate,  made  from  0.005  inch  (0.13  mm)  thick 
sheet  brass  shim  stock,  is  sandwiched  between  the  halves  of  the 
transducer  body  which  are  epoxied  and  screwed  together.  A  float- 
ing element   4  mm  square   is  cast  in  place  atop  the  cantilever 
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plate  using  an  epoxy  potting  compound.  The  square  shape  of  the 
floating  element  was  selected  to  minimize  the  degradation  in 
transducer  sensitivity  caused  by  sealing  the  gap  around  the 
element. 

The  transducer  is  then  mated  to  an  aluminum  adapter  machined 
to  accomodate  the  test  fixture  as  indicated  in  Figure  2.3.  An 
aluminum  cylinder  also  mates  with  the  adapter,  shielding  the 
transducer  and  its  electronics  from  electromagnetic  interference. 
Brass  shims  between  the  transducer  and  the  adapter  are  used  to 
position  the  floating  element  parallel  to  the  top  of  the  adapter 
with  the  top  of  the  floating  element  protruding  less  than  0.1  mm 
from  the  surface  of  the  adapter.  The  tops  of  the  adapter  and 
floating  element,  spray  painted  with  flat  white  enamel  before  as- 
sembly, were  lightly  sanded  with  fine  wet  sandpaper  to  reduce  the 
misalignment  of  the  floating  element  with  the  wall  to  9  ym.  By 
sanding  perpendicular  to  the  transducer's  sensitive  direction, 
motion  of  the  floating  element  and  curvature  of  its  upper  surface 
were  minimized.  The  gap  around  the  element  is  not  symmetric, 
ranging  from  0.1  mm  to  0.25  mm  with  a  mean  of  0.18  mm.  An  accel- 
erometer  mount  is  provided  on  the  side  of  the  transducer  assembly 
inside  the  cylinder  in  order  to  monitor  the  acceleration  of  the 
transducer  in  its  sensitive  direction. 

Table  II  summarizes  the  key  lengths  of  the  floating  element 
in  both  dimensional  form  and  in  inner  coordinates,  nondimension- 
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alized  on  the  viscous  length  A*  characteristic  of   the  boundary 
layer  in  which  the  transducer  was  ultimately  tested. 


FLOATING 
ELEMENT 


GAP 


TRANSDUCER 
BODY 

TRANSDUCER' 

MOUNTED 
ELECTRONICS 


TEST 
FIXTURE 
ADAPTER 


1  cm3 


ALUMINUM 
CYLINDER 

ELECTRICAL 
CONNECTIONS 


FIGURE  2.3    Transducer  -  Pictorial 


TABLE  II 
Floating  Element  Dimensions 


Dimension 

S 

+ 
5 

max. 

0. 25  mm 

18 

gap  min. 

0. 10  mm 

7 

mean 

0. 18  mm 

13 

Protrusion 

9  urn 

0.69 

Element 
Size 

4  mm 

290 



C   is  based  on  the  experimental  result 
U  =  1.3  m/s  at  U  =  30.5  m/s. 
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2.1.2  Transducer  Electronics 

The  electronics  associated  with  the  transducer  fall  into  two 
categories:  electronics  physically  mounted  beneath  the  phenolic 
transducer  assembly,  and  signal  conditioning  electronics  in- 
stalled in  a  separate  enclosure.  Each  set  of  electronics  uses 
shielded  precision  power  supplies  (Analog  Devices  model  904) 
which  are  bypassed  to  ground  through  1.0  PF  electrolytic  capaci- 
tors. Further,  each  integrated  circuit  power  supply  pin  is  by- 
passed to  ground  through  a  0.1  yF  ceramic  capacitor. 

Figure  2.4  highlights  the  significant  features  of  the  elec- 
tronic circuitry. 


R 


• — <vw* 


Floating 
Element 
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FIGURE  2.4   Circuitry  -  Block  Diagram 
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Transducer-mounted  electronics  consist  of  the  frequency 
modulator  and  a  buffer  amplifier.  The  frequency  modulator  is  an 
astable  multivibrator  which  produces  a  square  wave  nominally 
6.2  V  in  amplitude  at  a  frequency  determined  by  the  variable  ca- 
pacitor formed  by  the  fixed  and  cantilever  plates.  When  the  can- 
tilever plate  is  undeflected,  the  capacitance  measures  about 
30  pF  and  the  oscillator's  frequency  is  about  55.6  kHz.  An  LM310 
voltage  follower  is  used  as  an  isolation  amplifier  that  drives  a 
shielded  BNC  transmission  line  to  the  signal  conditioning  elec- 
tronics. Power  is  supplied  to  the  transducer-mounted  electronics 
through  a  shielded  cable  from  a  power  supply  in  a  separate  enclo- 
sure.  The  multivibrator  circuit  is  taken  from  Jung  (1980). 

The  signal  conditioning  electronics  are  mounted  in  a  shield- 
ed aluminum  box  containing  an  integral  power  supply.  An  Analog 
Devices  453L  precision  frequency-to-voltage  converter  demodulates 
the  transducer's  frequency  modulated  signal,  converting  it  to  a 
voltage  that  is  proportional  to  the  received  frequency.  Both  the 
gain  and  the  offset  of  the  frequency-to-voltage  conversion  are 
adjustable  for  calibration.  An  AD524  instrumentation  amplifier 
(with  selectable  gains  of  1,  10,  100,  or  1000)  amplifies  the  de- 
modulated signal  and  drives  both  an  active  filter  and  a  transmis- 
sion line  for  measuring  the  fluctuating  signal.  The  active 
filter,  taken  from  Lancaster  (1982),  is  a  first-order  low  pass 
filter  with  cutoff  frequency  of  0.2  Hz.  The  filter  drives  a  sep- 
arate line  for  measuring  the  steady  transducer  output. 


2.2  TRANSDUCER  RESPONSE  CHARACTERISTICS 


It  is  shown  in  Appendix  A. 3  that  the  response  of  the  trans- 
ducer is  described  by  the  following  equation  in  the  limit  of 
small  deflections  of  the  floating  element: 


3  A  G  f    _ 

V(t)   =  V   +  - —  t  (t)  (2.1) 

0      8  k  d 
s  o 


where 


x  (t)    =  the  mean  wall  shear  stress  acting  on  the   floating 

w 

element; 

V(t)     =  the  output  voltage  from  the   signal  conditioning 
electronics; 

A       =  the  surface  area  of  the  floating  element   exposed 
to  the  flow; 

d       ■  the  mean  separation  between  the  undeflected  capa- 
citor plates; 
f       =  the  oscillator  center  frequency; 

G       =  the  total  gain  in  the   frequency-to-voltage  con- 
verter and  the  instrumentation  amplifier; 

k      =  the  transverse  stiffness  of   the  cantilever  beam 
s 

when  statically  loaded  at  its  free  end;   and 
V      =  the  offset  voltage. 


Equation  2.1  is  applicable  for  small  displacements  of  the  float- 


^o 


ing  element  and  for  frequencies  much  smaller  than  the  first  reso- 
nant frequency  of  the  beam  element. 

Equation  2.1  may  be  rewritten  as 


V(t)   =  V   +  s   t  (t)   ,  (2.2) 

O        T   W 


where 


au      3  A  G  f 
-   o  v         e    o 

st  =  a7~  =  =  °-40  v/pa  (2.3) 

w      8  k  d 
s  o 


is  the  transducer's  sensitivity.   A  related   sensitivity   is   the 
transducer  output  per  micron  of  floating  element  deflection 

ks        3  G  f 
s    =  __§       o 

E     A  °x  =   0-56  V/ym   .  (?    ,. 

e        8  d  Kl.t*) 

o 


The  limits  of  the  validity  of  equation  2.1  may  be  investi- 
gated by  using  the  predicted  beam  element  deflection  shape  deter- 
mined in  Appendix  A. 3  to  evaluate  the  theoretical  transducer  out- 
put.  From  the  appendix, 


V(t)  -  Vo  1  (2.5) 

-   1  (A1.15) 


G  f 
o 


J  1  +  n(C,t) 
o 


and 


n(C,t)  = 


2d 


(U2     -  53)  , 


(2.6) 
(A3. 4) 


V(t)  -  V 

where is  the  normalized  transducer  output  andn(£,t)  is 

G  r 
o 

the  deflection  shape  of  the  beam  element  normalized  on  d 
Equation  2.6  describes  the  static  deflection  shape  of  the  beam 
element  for  a  force  applied  to  the  floating  element  and  is  valid 
for  frequencies  much  smaller  than  the  fundamental  mechanical  re- 
sonance. The  integral  in  equation  2.5  may  be  performed  numeri- 
cally using  equation  2.6  and  the  normalized  transducer  output 
plotted  as  a  function  of  the  normalized  end  deflection  of  the 
beam  -j —  .   As  shown  in  Figure  2.5,  the  transducer  is  linear  to 


within  1.25  percent  proportional  error  up  to 

2 

that  k,   =22.5  N/m,  A   =  16  mm  ,  and  d   - 
s  '     *     e  '      o 


=  0.05.   Given 


375  ym,  the  device 


is  effectively  linear  up  to 


w 


-  26  Pa  since 


w 


k 

k    d 

s 

i     i 

s     o 

F. 

sss 

■ 

£ 

=s 

A 

max 

A 

d 

max 

e 

e 

( 

max 


(2.7) 


For  the  turbulent  boundary  layers  considered  in  the  current   re- 


search, 


w 


max 


5  Pa  so  that  the  normal  operating  regime  is  the 


much  smaller  box  shown  in  Figure  2.5. 


The  natural  frequencies  of  the  transducer  response  are 
obtained  in  Appendix  B  by  evaluating  the  eigenvalues  of  the  beam 
element's  free  vibration  response.  For  the  first  two  vibration 
modes  of  the  beam  element,  the  rotational  inertia  of  the  floating 
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element  was  calculated  to  affect  the  natural  frequencies  by  less 
than  1.0  percent  and  so  was  neglected.  Calculation  of  the  trans- 
ducer's natural  frequencies  is  complicated  by  an  ambiguity  in  the 
beam  length  introduced  by  the  spatial  extent  of  the  floating  ele- 
ment. Contrary  to  the  model,  the  floating  element  is  not  a  point 
mass  but  rather  extends  over  6  percent  of  the  cantilever  plate's 
length.  The  fundamental  resonant  frequency  of  the  beam  element 
is  estimated  from  the  model  to  be  61.7  Hz,  with  an  uncertainty  of 
6.3  percent.  Measurements  of  the  dynamic  characteristics  of  the 
transducer  presented  in  greater  detail  in  section  4.1  indicate 
that  the  fundamental  resonance  is  actually  63  Hz. 

Characteristics  of  the  transducer  are   summarized   in  Table 
III. 


TABLE  III 
Transducer  Characteristics 


Mass:   Floating  Element      76  mg 
Cantilever  Plate     480  mg 

=  16  g/m 

Stiffness:  22.5  N/m 

Fundamental  Resonance:       63  Hz 
Sensitivity:  s  =0.40  V/Pa 

T 


s  =  0. 56  V/ym  end 

deflection 


Spatial  Resolution:        4  mm  x  4  mm 


JO 


Inclination  of  the  transducer  from  the  vertical  in  a  gravity 
field  generates  a  response  from  the  transducer,  since  then  a  com- 
ponent of  the  weight  of  the  beam  element  is  transverse  to  the 
beam.  The  transducer  exhibits  buckling  behavior  when  inclined, 
because  the  transverse  force  tends  to  increase  the  deflection  of 
the  beam  and  itself  increases  as  the  local  inclination  angle  of 
the  beam  increases.  Transducer  inclination  response  is  analyzed 
in  Appendix  A. 4  in  the  limit  of  small  static  inclinations  of  the 
transducer  base.  The  expression  obtained  in  the  appendix  for  the 
deflection  of   the  beam  element  n(£,t)with  inclination  angle 

0  (equation  A4.7)  is  then  inserted  in  equation  2.5,  which  is  in 

o 

turn  integrated  numerically  to  obtain  the  transducer  output  as  a 
function  of  inclination  angle  plotted  in  Figure  2.6. 

The  theoretical  response  of  the  transducer  to  inclination  is 
sensitive  to  the  beam  length  L  and  mean  capacitor  spacing  d  ,  as 
indicated  in  Figure  2.6  by  comparing  the  solid  and  dashed  lines. 
The  solid  line  is  the  theoretical  transducer  response  for  the 
beam  length  and  mean  capacitor  spacing  used  in  all  calculations; 
the  dashed  line   is   the  theoretical  response  for  a  beam  length 

1  mm  (3.2%)  smaller  and  mean  capacitor  spacing  25  ym  (6.7%) 
larger.  The  variations  in  the  parameters  used  to  generate  the 
dashed  line  are  representative  of  uncertainties  in  those  parame- 
ters due  to  manufacturing  limitations. 

A  simple  experiment  was  conducted  to  verify  the   theoretical 


Jl 


response  predicted  in  Figure  2.6.  Experimental  results  are  also 
plotted  in  Figure  2.6  and  may  be  seen  to  lie  within  the  limits  of 
uncertainty  in  predicted  response.  Also  indicated  on  the  figure 
are  the  expected  operating  limits  of  the  transducer  in  the  cur- 


T 
W 


=  5.0  Pa  at 


rent  research  from  Figure  2.5,  corresponding  to 
zero  inclination  angle.  The  transducer  output  due  solely  to  ro- 
tation of  the  transducer  base  one  degree  from  the  vertical  is  of 
the  same  order  as  the  transducer  output  due  solely  to  wall  shear 
stress,  indicating  the  importance  of  maintaining  the  attitude  of 
the  transducer. 
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3.  EXPERIMENTAL  EQUIPMENT  AND  INSTRUMENTATION 


3.1  EXPERIMENTAL  FACILITY 

The  transducer  was  tested  in  the  MIT  low-noise, 
low-turbulence  open  circuit  wind  tunnel.  This  facility  is  shown 
in  Figure  3.1  and  is  described  in  detail  by  Hanson  (1969).  It 
has  a  flow  speed  range  of  20  to  50  m/s,  controlled  by  varying  the 
speed  of  a  DC  motor  driving  the  blower.  The  array  of  turbulence 
reducing  screens  between  the  honeycomb  flow  straightener  and  the 
settling  zone  maintains  the  free  stream  turbulence  level  less 
than  approximately  0.05%.  A  test  section  38  cm  (15  inches) 
square  and  183  cm  (6  feet)  long  was  mounted  in  an  open  jet  confi- 
guration within  the  airtight  blockhouse  immediately  behind  the 
contraction.  The  test  section,  lined  with  smooth  sound  absorp- 
tive material,  has  movable  walls  for  controlling  the  longitudinal 
pressure  gradient.  A  mechanical  x-y  traverse  accurate  to  0.025 
mm  (0.001  inch)  vertically  and  to  2.5  mm  (0.1  inch)  horizontally 
is  installed  on  the  test  section  and  used  to  position  a  hot  wire 
or  a  pressure  probe  during  the  measurement  of  mean  flow  proper- 
ties in  the  test  section. 
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3.2  TEST  FIXTURE 

The  test  fixture  is  a  damped  aluminium  flat  plate  described 
in  detail  by  Gedney  (1979),  who  used  it  while  investigating  tran- 
sition region  wall  pressure  fluctuations.  The  plate,  sketched  in 
Figure  3.2,  is  38  cm  (15  inches)  wide  and  122  cm  (48  inches) 
long;  the  upper  surface  is  polished  to  a  0.1  micron  root  mean 
square  (rms)  surface  roughness.  The  plate's  leading  edge  is  a 
solid  tool  and  jig  plate,  1.27  cm  (0.50  inch)  thick,  machined 
into  a  6:1  ellipse.  The  test  plate's  trailing  edge  consists  of  a 
thin  splitter  plate  which  extends  rearward  to  prevent  coherent 
vortex  shedding  from  creating  plate  vibrations. 

The  test  plate  was  mounted  horizontally  in  the  test  section 
13  cm  (5  inches)  above  the  bottom  to  minimize  secondary  flow  ef- 
fects observed  by  Shapiro  (1977).  The  walls  of  the  test  section 
were  then  adjusted  to  provide  a  minimum  pressure  gradient.  The 
variation  in  pressure  obtained  over  the  length  of  the  plate  at  a 
free  stream  velocity  of  35  m/s  is  shown  in  Figure  3.3.  The  maxi- 
mum variation  observed  (excluding  the  leading  edge)  was  less  than 
2%  of  the  total  head.  The  transducer  was  installed  at  the  loca- 
tion indicated. 
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FIGURE  3.3  Longitudinal  Pressure  Variation 
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A  flow  visualization  to  check  the  effect  of  transverse  wall 
turbulence  was  performed  using  Shapiro's  technique,  also  used  by 
Gedney.  The  results  were  identical  to  Gedney's  as  shown  in  Fig- 
ure 3.4.  The  shear  stress  transducer  was  tested  in  the  7.6  cm  (3 
inch)  hole  at  the  location  indicated. 

The  boundary  layer  was  tripped  to  initiate  turbulence  by  two 
thicknesses  of  duct  tape  1  cm  in  width  and  0.5  mm  in  total 
height.  The  tape  was  placed  on  the  upper  surface  of  the  plate 
2  cm  from  and  parallel  to  the  leading  edge.  According  to  the 
criterion  developed  by  Gibbings  [White  (1976)],  fully  effective 
tripping  to  turbulent  flow  is  assured  if 


U  k 

-V  "   826  (3-D 


where  k  is  the  height  of  the  surface  roughness.  In  the  present 
work,  k  -  0.5  mm  ensures  immediate  transition  for  free  stream 
velocity  greater  than  25  m/s. 
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3.3  FLOATING  ELEMENT  MISALIGNMENT 


Allen  (1977)  has  indicated  that  the  major  sources  of  error 
in  floating  element  shear  stress  measuring  devices  are  lip  forces 
and  normal  forces  which  contaminate  the  friction  force  that  the 
device  is  being  used  to  measure.  Lip  forces  are  developed  by  a 
nonuniform  pressure  distribution  around  the  periphery  of  the 
floating  element  caused  by  either  flow  through  the  gap  around  the 
element  or  misalignment  of  the  floating  element  with  respect  to 
the  wall.  Normal  forces  introduce  error  when  a  nonuniform  pres- 
sure distribution  over  the  surface  of  the  floating  element  shifts 
the  normal  force's  center  of  effort,  generating  a  moment  which 
deflects  the  floating  element. 

Allen  found  that  errors  in  the  measured  wall  shear  stress 
due  to  lip  and  normal  forces  can  be  larger  than  the  actual  wall 
shear  stress  by  as  much  as  a  factor  of  4.  His  results  indicate 
that  proper  alignment  of  the  floating  element  in  the  wall  is  the 
key  factor  in  minimizing  error,  a  finding  corroborated  by  Burke 
(1984)  who  used  a  floating  element  to  measure  skin  friction  in 
the  MIT  water  tunnel. 

Proper  alignment  of  the  floating  element  in  the  wall  is 
therefore  critical  in  the  current  design,  and  was  assured  by 
measurement  with  a  KD-45  Fotonic  Sensor  made  by  Mechanical  Tech- 
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nology  Incorporated.  The  Fotonic  Sensor  is  a  non-contact  measur- 
ing device  that  projects  several  beams  of  light  through  fiberop- 
tic strands  and  then  measures  the  quantity  of  light  reflected 
from  an  object  to  determine  the  distance  to  the  object.  The  Fo- 
tonic Sensor  was  used  in  its  less  sensitive  range  (  ±2  ym)  to  de- 
termine that  the  floating  element  protruded  9  ym  from  the  wall 
(y   =  0.69),  which  was  judged  to  be  satisfactory. 


3.4  INCLINATION  RESPONSE 


Experimental  verification  of  the  transducer's  response  to 
inclinations  from  the  vertical  were  made,  as  shown  in  Figure  2.6. 
The  transducer  was  attached  to  the  mid-span  of  a  sturdy  1.83  m 
long  aluminum  I-beam  which  was  raised  or  lowered  at  one  end  a 
known  distance  to  generate  inclination  angles  smaller  than  one 
degree.  An  inclination  angle  of  one  degree  is  obtained  when  one 
end  of  the  I-beam  is  32  mm  higher  than  the  other  end,  so  it  is 
possible  to  maintain  accuracy  to  about  0.01  degree.  The  effect 
of  I-beam  flexure  on  inclination  angle  is  small  to  begin  with  and 
further  reduced  by  placing  the  transducer  at  the  I-beam's 
mid-span,  so  that  I-beam  flexural  effects  have  been  neglected. 
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3.5  INSTRUMENTATION 


3.5.1  Dynamic  Characteristics 

Transducer  dynamic  characteristics  were  investigated  using 
the  instrumentation  shown  in  Figure  3.5.  The  transducer  was 
mounted  in  a  phenolic  bracket  and  secured  to  the  vibration  ex- 
citer so  that  the  transducer  could  be  shaken  either  in  the  plane 
of  the  wall  or  normal  to  the  wall.  As  shown  in  Appendix  A. 2,  the 
transducer's  response  to  acceleration  of  its  base  and  its  re- 
sponse to  a  shear  force  acting  on  the  floating  element  are  dynam- 
ically equivalent,  although  caution  must  be  exercised  in  inter- 
preting results  since  the  acceleration  of  the  transducer's  base 
also  introduces  an  effective  force  along  the  length  of  the  canti- 
lever plate  which  is  absent  during  normal  operation  of  the  trans- 
ducer. An  accelerometer ,  calibrated  using  a  General  Radio  model 
1557-A  vibration  calibrator,  was  attached  to  either  the  phenolic 
bracket  or  to  the  transducer  body  to  measure  the  vibration  of  the 
transducer  body.  A  Hewlett-Packard  5423A  structural  dynamics  an- 
alyzer was  used  to  generate  the  random  noise  excitation  of  the 
vibration  exciter  and  also  to  calculate  the  transfer  function 
between  the  transducer  output  and  transducer  body  acceleration. 
The  spectrum  of  the  random  noise  excitation  was  shaped  using  an 
MXR  31-band  equalizer  so  that  the  spectrum  of  the  transducer  body 
acceleration  was  approximately  uniform  over  the  frequency  band  of 
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interest 


3.5.2  Test  Conditions 

The  free  stream  velocity  into  the  test  section  was  monitored 
by  measuring  the  pressure  drop  across  the  wind  tunnel  contraction 
section  with  a  Betz  micromanometer .  A  calibration  was  performed 
by  Shapiro  (1977)  to  ensure  an  accuracy  of  0.1  m/s. 

The  longitudinal  pressure  variation  along  the  test  plate  was 
measured  using  a  small  static  pressure  tube  (United  Sensor  and 
Control  Corp.  PSA-16)  and  a  Validyne  DP15-24  differential  pres- 
sure transducer.  A  Validyne  model  CD23-1071  electronic  manometer 
powered  the  pressure  transducer  and  provided  a  digital  readout  of 
the  differential  pressure  in  millimetres  of  water. 

Mean  velocity  profile  measurements  were  performed  using  a 
hot  wire  anemometer  probe  ( The rmo-Sys terns  Inc.  probe  model 
1261-T1.5)  with  a  Disa  model  55005  battery-powered  constant  tem- 
perature anemometer  and  a  linearizer  built  by  Moeller  using  an 
Analog  Devices  433B  programmable  function  generator.  The  hot 
wire  probe  was  calibrated  against  a  pitot  tube  used  with  the  Val- 
idyne DP15-24  differential  pressure  cell.  Mean  velocity  was 
measured  on  a  digital  voltmeter  after  being  filtered  by  a  0.1  Hz 
active  low  pass  filter  built  by  Lueptow  (1984). 
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3.5.3  Transducer  Calibration 

Wall  shear  stress  transducer  electronics  were  calibrated 
using  a  General  Radio  model  1210-C  function  generator  in  conjunc- 
tion with  a  Hewlett-Packard  model  5532A  electronic  counter  and  a 
digital  voltmeter.  With  gain  of  the  AD524  instrumentation  am- 
plifier set  at  10,  the  gain  of  the  frequency  to  voltage  converter 
was  adjusted  so  that  a  change  in  input  frequency  of  500  Hz  (from 
55.6  kHz  to  either  55.1  kHz  or  56.1  kHz)  produced  a  change  in 
output  voltage  of  approximately  5.0  V.  The  exceptional  sensitiv- 
ity of  the  transducer  made  workbench  calibration  of  the  entire 
transducer  unfeasible,  so  the  transducer  was  calibrated  against 
the  law  of  the  wall  in  situ. 


3.5.4  Wall  Shear  Stress 

The  wall  shear  stress  was  measured  using  the  instrumentation 
of  Figure  3.6.  Coherence  between  the  acceleration  of  the  trans- 
ducer and  transducer  output  was  monitored  using  the  structural 
dynamics  analyzer  to  ensure  that  transducer  acceleration  did  not 
significantly  affect  transducer  output.  The  structural  dynamics 
analyzer  was  also  used  in  its  time  record  mode  to  make 
long-period  zero  drift  measurements. 
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EXPERIMENTAL  RESULTS  AND  DISCUSSION 


4.1  DYNAMIC  CHARACTERISTICS 

Transducer  dynamic  characteristics  were  evaluated  for  two 
different  transducer  configurations:  one  with  a  seal  in  the  gap 
surrounding  the  floating  element  and  another  with  the  gap 
unsealed.  For  brevity,  these  two  configurations  will  be  referred 
to  as  the  "sealed"  and  "unsealed"  configurations,  respectively. 
Investigation  of  the  dynamic  characteristics  of  the  two  gap  con- 
figurations was  conducted  in  two  phases.  In  the  first  phase,  the 
accelerometer  used  to  measure  transducer  body  acceleration  was 
attached  to  the  phenolic  bracket  mounting  the  transducer  to  the 
vibration  exciter.  This  test  configuration  permitted  the  rota- 
tion of  the  transducer  in  the  bracket  and  was  used  to  verify  that 
the  transducer  can  discriminate  against  signals  orthogonal  to  its 
sensitive  axis  in  the  plane  of  the  wall,  but  introduced  errors  in 
that  the  accelerometer  was  not  mounted  directly  on  the  transducer 
body.  The  second  phase  utilized  an  accelerometer  attached  di- 
rectly to  the  transducer  body  in  the  location  indicated  on  Figure 
2.2.  In  this  later  phase,  attention  was  devoted  exclusively  to 
improving  understanding  of  the  dynamics  of  the  transducer  in  its 
sensitive  direction,  the  orthogonality  of  the  device  having 
already  been  demonstrated. 
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Phase  1.  The  sealed  configuration  in  the  initial  phase  used 
as  the  gap  sealing  material  an  epoxy,  ECCOGEL  1365-15,  made  by 
Emerson  and  Cuming.  Measuring  only  15  on  the  Shore  hardness 
scale  commonly  used  in  the  polymer  industry,  the  epoxy  is  very 
soft  but  is  nearly  incompressible.  Owing  to  the  manner  of  the 
transducer's  construction,  it  was  not  possible  to  determine  the 
geometry  of  the  gap  seal  beneath  the  surfaces  of  the  adapter  and 
floating  element  after  the  epoxy  cured. 

The  transfer  function  between  the  transducer  output  and  the 
transverse  acceleration  of  the  transducer  body  is  shown  in  Figure 
4.1  (unsealed  configuration)  and  Figure  4.2  (sealed  configura- 
tion). The  magnitude  plots  are  normalized  on  the  peak  magnitude 
of  the  unsealed  transducer's  transfer  function  in  order  to  high- 
light the  order  of  magnitude  loss  in  response  caused  by  the  epoxy 
seal  in  the  gap.  The  fundamental  resonant  frequency  of  the 
unsealed  transducer  (Figure  4.1)  is  approximately  63  Hz,  very 
close  to  the  61.7  Hz  predicted  during  design.  Figure  4.2  indi- 
cates that  the  addition  of  the  epoxy  seal  leads  to  an  increase  in 
the  fundamental  resonant  frequency  to  about  360  Hz  as  well  as  an 
apparent  reduction  in  sensitivity. 

The  apparent  loss  of  sensitivity  is  misleading.  Table  V 
(Appendix  B)  lists  the  resonant  frequencies  for  the  first  three 
vibration  modes  of  the  beam  element  subjected  to  clamped-free, 
clamped-pinned,  and  clamped-clamped  boundary  conditions.   The  ob- 
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served  fundamental  resonant  frequency  of  the  sealed  transducer 
(360  Hz)  is  bracketed  by  that  of  a  clamped-pinned  beam  element 
(341  Hz)  and  a  clamped-clamped  beam  element  (494  Hz),  neglecting 
the  rotational  inertia  of  the  floating  element  which  was  found  to 
be  insignificant.  This  suggests  that  the  beam  element  in  the 
sealed  configuration  exhibits  a  pinned  condition  at  the  floating 
element  or,  more  correctly,  a  mixed  condition  that  is  a  hybrid  of 
pinned  and  clamped  conditions  but  closer  to  pinned.  If  such  is 
the  case,  then  the  sealed  transducer  does  not  respond  to  wall 
shear  stress,  since  the  floating  element  is  located  at  a  node  of 
the  beam  element's  response  and  therefore  never  translates.  The 
sealed  transducer  is  thus  incapable  of  responding  to  forces  ap- 
plied to  the  floating  element,  but  rather  responds  only  to  forces 
applied   to   the  cantilever  plate.   This  is  shown  in  Figure  4.2, 

where  the  transducer  is  responding  to  the  D'Alembert  or   inertia 

d2z 

force  _   u  °  (using  the  notation  of  Appendix  A. 2)  produced  on 

dt2 

the  cantilever  plate  by  the  transverse  acceleration  of  the  trans- 
ducer base.  The  use  of  a  seal  in  the  gap  around  the  floating 
element  is  self-defeating  in  this  instance,  eliminating  transduc- 
er response  to  the  wall  shear  stress. 

The  resonances  of  both  the  sealed  and  unsealed  configura- 
tions exhibit  considerable  damping,  approximately  14  percent  of 
critical  damping  in  each  case.  The  unsealed  configuration's 
response  in  Figure  4.1a  shows  an  irregularity  at  60  Hz  caused  by 
noise  from  the  field  windings  of  the  vibration  exciter,  while  the 


JJ 


poor  response  below  10  Hz  is  due  to  a  low  frequency  limitation  of 
the  vibration  exciter.  The  phase  angle  of  both  configurations' 
transfer  function  is  anomalous,  since  it  does  not  change  by  180 
degrees  as  it  should  when  passing  through  an  isolated  resonance. 
The  irregularity  in  Figure  4.1a  near  50  Hz  and  the  similar  irre- 
gularity in  Figure  4.2a  near  325  Hz  are  also  strange,  but  may  be 
due  to  the  vibration  of  the  bracket  securing  the  transducer  to 
the  exciter. 


Phase  2.  The  second  phase  of  the  investigation  of  transduc- 
er dynamics  was  undertaken  to  resolve  some  of  the  issues  raised 
concerning  the  phase  of  the  transfer  function  and  the  role  of  the 
mounting  bracket  in  introducing  the  anomalies  encountered  in  the 
earlier  phase.  Also  of  interest  was  confirmation  of  the  hypo- 
thesis that  the  floating  element  in  the  sealed  configuration  was 
effectively  held  motionless  by  the  seal.  In  this  phase,  the  ac- 
celerometer  used  to  measure  transducer  transverse  acceleration 
was  mounted  directly  on  the  phenolic  transducer  assembly  rather 
than  on  the  bracket  mounting  the  transducer  to  the  vibration  ex- 
citer. For  the  "sealed"  configuration,  a  single  piece  of  stan- 
dard vinyl  electrical  tape  was  securely  attached  to  the  surfaces 
of  the  floating  element  and  adapter  to  approximate  a  pinned  boun- 
dary condition  at  the  floating  element. 

Figures   4.3  and  4.4   illustrate   the   transfer   function 
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obtained  between  the  transducer  output  and  the  transverse  accel- 
eration of  the  transducer  body.  The  units  of  the  ordinate  in  the 
magnitude  plot  are  volts  per  g,  where  1  g  =  9.8  m/s  is  the  ac- 
celeration due  to  gravity. 

The  unsealed  configuration's  transfer  function  (Figure  4.3) 
is  virtually  identical  to  Figure  4.1,  with  the  following  impor- 
tant exceptions: 

-  the  phase  is  well-behaved  and  very  properly  shifts 
through  180  degrees  when  passing  through  the  resonance; 

-  the  anomaly  in  the  magnitude  near  50  Hz  has  vanished, 
as  expected;   and 

-  the  irregularity  in  magnitude  at  60  Hz  is  sharper  than 
that  encountered  earlier,  and  gives  rise  to  slight  ir- 
regularities at  120  Hz  and  180  Hz. 

As  discussed  earlier,  the  noise  at  60  Hz  (and  multiples  thereof) 
is  due  to  the  field  windings  in  the  vibration  exciter.  The 
slight  irregularity  in  magnitude  at  100  Hz  is  due  to  noise  from 
the  vibration  exciter's  blower  motor. 

The  "sealed"  configuration's  transfer  function  is  remarkably 
similar  to  that  of  Figure  4.2,  despite  the  difference  between  the 
two  seals  (electrical  tape  as  opposed  to  epoxy) .   In  Figure  4.4: 

-  the  fundamental  resonant  frequency  has  increased  to 
375  Hz  (from  360  Hz  in  Figure  4.2); 

-  the  damping  has  decreased  to  7  percent  of  critical   (in 
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contrast  to  14  percent  in  Figure  4.2); 

-  the  large  irregularity  in  magnitude   just   below  reso- 
nance has  vanished;   and 

-  the  phase  is  well-behaved  and  shifts   through   180  de- 
grees through  the  resonance. 

The  results  of  this  phase  of  the  investigation  indicate  that 
the  anomalies  encountered  were  indeed  experimental  and  confirm 
the  hypothesis  that  the  gap  seal  immobilizes  the  floating  element 
in  translation. 


4.2  WALL  SHEAR  STRESS 

4.2.1  Law  of  the  Wall 

The  law  of  the  wall  was  used  to  obtain  estimates  of  the  mean 
wall  shear  stress  for  comparison  with  the  steady  mean  wall  shear 
stress  measured  by  the  transducer.  Mean  velocity  profiles  were 
measured  above  the  floating  element  with  a  hot  wire  probe  at  free 
stream  velocities  of  26.0  m/s,  30.5  m/s,  and  34.5  m/s  correspond- 
ing to  momentum  thickness  Reynolds  numbers  of  3140,  3500,  and 
3760.  The  velocity  profiles  were  analyzed  using  the  FORTRAN  pro- 
gram BLPR0F2  in  which  the  user  iteratively  guesses  the  friction 
velocity  U  and  the  location  of  the  wall  (in  the  measured  coor- 
dinates) until  the  measured  data  coincides  in  the  inner  region  of 
the  profile  with  Musker's  explicit  formula  for  the  turbulent  mean 
velocity  profile  [see  Musker  (1979)].  BLPR0F2  permits  the  user 
to  fair  the  measured  data  into  the  wall  using  Musker's  formula 
and  then  integrates  the  resulting  faired  profile  using  trapezoi- 
dal rule  to  obtain  displacement  thickness  6*  and  momentum  thick- 
ness 6  . 

The  three  mean  velocity  profiles  are  plotted  in  Figure  4.5 
nondimensionalized  on  inner  parameters  of  the  boundary  layer.  It 
can  be  seen  that  the  velocity  profile  for  Re  =  3500  contains 
little  or  no  wake,  whereas  the  other  two  profiles  exhibit  a  no- 
ticeable wake.   The  velocity  profile  at  Re  =  3500  was  measured 
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during  a  different   test   period   than  the  velocity  profiles  at 

Re  -  3140  and  3760,  highlighting  the  difficulty  of   reproducing 

0 

test  conditions  exactly.  Since  the  anomaly  in  the  wake  at 
Re9  =  3500  does  not  affect  the  logarithmic  overlap  region,  the 
change  in  the  outer  flow  causing  the  different  wake  is  expected 
to  have  little  effect  on  the  wall  shear  stress.  Because  all 
three  mean  velocity  profiles  collapse  to  the  classical  turbulent 
mean  velocity  profile,  it  is  assumed  that  each  boundary  layer  ex- 
hibits fully  developed  turbulence  despite  being  tripped  into  tur- 
bulence so  that  the  application  of  the  law  of  the  wall  is  valid. 


4.2.2  Transducer  Measurement 

Steady.  Steady  mean  wall  shear  stress  was  measured  by  the 
transducer  at  the  same  Reynolds  numbers  as  the  velocity  profiles; 
the  results  are  listed  in  Table  IV,  where  they  are  compared  to 
the  results  obtained  from  the  law  of  the  wall.  Both  the  law  of 
the  wall  data  and  the  transducer-measured  data  are  plotted  in 
Figure  4.6,  taken  from  skin  friction  data  analyzed  by  Murlis, 
et.  al.  ,  (1982).  Transducer  measurements  at  Re0  -  3140  and  3500 
are  excellent,  while  that  for  Re  -  =  376Q  is  a  bit  high  but  still 
within  experimental  uncertainty,  estimated  in  Appendix  C  to  be  66 
percent.  The  acceleration  of  the  transducer  was  neglected  in  an- 
alyzing transducer  output,  since  low  coherence  between  the  two 
(at  most   10%  to  20%)  indicated  that  there  was  little  systematic 
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relationship  between  the  transducer  output  and  acceleration  of 
the  transducer  body. 


TABLE  IV 
Results 


Ree 

3140 

350C 

3760 

x 
law  of  the  wall 

1.25  Pa 

2.05  Pa 

2.20  Pa 

measured 

1.3  Pa 

1.7  Pa 

3.1  Pa 

Transducer  zero  drift  was  surprisingly  low  despite  the 
transducer's  extreme  sensitivity,  averaging  20%  of  the  mean  sig- 
nal over  the  period  of  the  measurement  (1  to  1.5  minutes).  Zero 
drift  was  measured  by  covering  the  floating  element  with  an  index 
card  while  the  wind  tunnel  was  running  at  the  proper  speed.  A 
short  length  of  wire  0.5  mm  in  diameter  was  taped  beneath  the 
card  as  a  spacer  to  keep  the  card  from  contacting  the  floating 
element.  Sample  plots  of  the  zero  drift  measured  by  the  struc- 
tural dynamics  analyzer  over  a  period  of  1  minute  are  shown  in 
Figure  4.7  where  they  are  compared  to  the  mean  output  voltage  de- 
termined by  subtracting  the  average  reference  voltage  from  the 
transducer  total  output  voltage.  Each  data  run  took  3  to  4 
minutes  from  the  start  of  the  initial  zero  measurement  to  the 
measurement  of  the  wall  shear  stress  and  the  concluding  zero 
measurement . 
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FIGURE  4.6. 


Skin  friction  coefficient  variation  with  momentum  thick- 
ness Reynolds  number  [from  Murlis,  et.  al. ,  (1982)1  :  +, 
measured  data  using  wall  shear  transducer;  o,  data  from 
law  of  the  wall. 
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FIGURE  4.  7   Zero  Drift 
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Unsteady.  The  unsteady  mean  wall  shear  stress  spectrum  for 
Reg  =  3140  and  3760  is  shown  with  the  background  spectrum  in  Fig- 
ure 4.8.  Although  the  plots  show  some  dynamic  character  of  the 
unsteady  wall  shear  stress,  they  show  little  of  its  structure 
owing  to  the  low  resonant  frequency  of  the  transducer.  If  it  is 
assumed  that  an  instrument  is  usable  up  to  a  frequency  one  order 
of  magnitude  less  than  its  first  natural  frequency,  then  the 
transducer  can  only  be  used  up  to  6.3  Hz.  Burton  (1974)  indi- 
cates that  the  frequency  related  to  sublayer  bursting  in  a  turbu- 
lent boundary  layer  is  given  approximately  by 


^   ■  0.03  (4.D 

00 


where  6*  is  the  displacement  thickness.  For  the  conditions  under 
which  the  transducer  was  tested,  equation  4.1  predicts  the  burst 
frequency  to  be  about  300  to  500  Hz,  well  above  the  nominal 
bandwidth  of  the  transducer  and  thus  unmeasurable. 
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FIGURE  4.8    Wall  Shear  Stress  Spectral  Density 
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5.  RECOMMENDATIONS 


In  the  short  term,  the  wall  shear  transducer  is  usable  for 
measuring  the  steady  mean  wall  shear  stress.  The  sensitivity  of 
the  device  is  so  great,  however,  that  calibration  must  be  per- 
formed in  situ,  preferably  in  a  laminar  boundary  layer  where  the 
wall  shear  stress  is  well  known.  The  electronics  should  be  up- 
graded by  the  installation  of  precision  instrumentation 
amplifiers  to  replace  the  commercial-grade  amplifiers  currently 
in  use.  This  should  reduce  the  high  frequency  noise  observed  at 
the  transducer  output  and  improve  the  stability  of  the  astable 
multivibrator,  whose  center  frequency  was  observed  to  drift  from 
55.6  kHz  to  46.5  kHz  over  seven  days  of  continuous  operation. 


Long  range  plans  for  the  transducer  should  include  increas- 
ing its  bandwidth,  decreasing  its  sensitivity  to  inclination,  and 
gaining  an  understanding  of  how  the  gap  around  the  floating  ele- 
ment affects  the  shear  stress  measurement  (given  that  the  gap 
cannot  be  sealed,  as  the  current  research  suggests). 

As  regards  the  bandwidth,  the  first  resonant  frequency  of 
the  beam  element  needs  to  be  increased  at  least  two  orders  of 
magnitude  so  that  the  transducer  can  measure  frequencies  near  the 
burst   rate.   A  larger  bandwidth  is  desirable  so  that  the  device 
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can  always  be  used  more  than  an  order  of  magnitude  below  its  fun- 
damental resonance,  which  is  far  less  complicated  than  attempting 
to  use  an  instrument  through  its  first  resonance.  The  transducer 
should  be  as  simple  to  use  as  possible.  The  linearity  of  the 
transducer,  coupled  with  sufficient  bandwidth,  will  make  it  a 
powerful  instrument  for  measuring  the  dynamics  of  the  wall  shear 
stress.  As  linear  system  theory  demonstrates,  only  in  a  linear 
system  can  the  statistics  of  an  output  (such  as  transducer  output 
voltage)  be  related  definitively  and  unambiguously  with  statis- 
tics of  the  input  (such  as  wall  shear  stress). 

Next,  the  sensitivity  of  the  transducer  to  inclinations  from 
the  vertical  should  be  reduced  by  at  least  two  orders  of  magni- 
tude so  that  inclination  has  a  second-order  effect  on  transducer 
output  (it  is  assumed  that  the  steady  transducer  signal  is 
zeroth-order  while  the  unsteady  signal  is  first-order) .  Coupled 
with  the  need  to  maintain  roughly  the  same  sensitivity  to  wall 
shear  stress  s  ,  the  requirements  for  an  improved  physical  design 
imply  that  the  mass  of  the  floating  element  and  capacitor  plate 
must  be  reduced  several  orders  of  magnitude.   If  the  size  of   the 

floating  element   is   reduced,   then  equation  2.3  indicates  that 

f   G 

must  be  increased  a  compensating  amount  to  maintain  the 


ks  do 


desired  s  .   Decreases  in  k  require  further  decreases  in  mass  in 

T  S 

order  to  maintain  the  required  bandwidth.  Increases  in  G  lead  to 
only  marginal  improvements  in  transducer  performance  since  noise 
is  generally  amplified  along  with  the  signal. 
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Finally,  and  most  important,  a  better  understanding  of  the 
role  of  the  floating  element  gap  in  affecting  the  boundary  layer 
must  be  reached.  There  is  little  in  the  literature  that  sheds 
any  light  on  this  issue.  Allen  (1977)  identified  that  small  gaps 
increase  the  sensitivity  of  the  floating  element  transducer  to 
errors  introduced  by  misalignment  of  the  floating  element  with 
the  wall.  Winter  (1977)  mentions  that  the  effect  of  gaps  perpen- 
dicular to  the  flow  is  poorly  understood,  but  that  the  effect  of 
the  gap  should  be  negligible  if  the  gap  width  is  smaller  than  100 
in  wall  coordinates.  McMichael,  et.  al. ,  (1979)  found  that  a 
transverse  slot  in  the  wall  of  a  turbulent  flow  affects  the  boun- 
dary layer  downstream  of  the  slot.  In  the  authors'  words,  "Evi- 
dently, temporary  release  of  the  no  slip  condition  as  the  bounda- 
ry layer  passes  over  the  slot  results  in  a  local  decrease  in  tur- 
bulence production  as  indicated  by  the  weakened  intensity  near 
the  wall  just  downstream  of  the  slot."  (p.  417).  Since  the  gap 
is  adjacent  to  the  wall  area  where  the  shear  stress  will  be  meas- 
ured, it  is  unclear  whether  or  not  the  boundary  layer  will  "heal" 
Itself  sufficiently  to  permit  an  unbiased  measure  of  the  wall 
shear  stress,  particularly  if  the  measurement  area  is  small.  The 
effect  of  the  gap  is  a  fundamental  issue  which  has  not  yet  been 
adequately  addressed. 


The  most  fruitful  course  of  action  appears  to  be  the  micro- 
miniaturization of   the  transducer  to  a  scale  where  a  very  small 
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gap  and  precise  alignment  of  the  floating  element  can  be  mainta- 
ined. Peterson  (1982)  discusses  the  merits  of  unicrystalline 
silicon  dioxide  at  great  length;  small  cantilever  beams  8  ym 
long,  2  ym  wide,  and  0.1  ym  thick  with  a  fundamental  resonant 
frequency  of  1.2  MHz  have  already  been  made  from  the  material.  A 
microminiaturized  wall  shear  stress  transducer  with  such  a  small 
cantilever,  using  the  capcitance-measuring  technique  outlined  in 
the  current  research,  has  a  number  of  distinct  advantages.  A 
floating  element  with  sensing  area  2  ym  by  0.1  ym  would  have  su- 
perb spatial  resolution,  excellent  sensitivity,  and  a  large 
bandwidth.  The  considerable  reduction  in  mass  would  greatly  re- 
duce the  sensitivity  of  the  device  to  inclination.  Further,  it 
should  be  possible  to  keep  the  gap  around  the  floating  element  on 
the  order  of  about  10  angstroms  so  that  the  gap  has  a  negligible 
effect  on  the  flow.  Precise  alignment  of  the  floating  element  is 
possible  on  the  microscopic  scale,  which  may  circumvent  the  ex- 
treme sensitivity  of  the  floating  element  to  misalignment-induced 
errors  when  the  gap  is  small.  Recent  trends  suggest  that  such  a 
transducer  could  be  mounted  on  a  20  ym  chip  with  its  own  acceler- 
ometer  and  compensation  network,  permitting  an  array  of  such 
transducers  to  measure  the  spatial  distribution  of  the  shear 
stress  over  the  surface  of  the  wall. 
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6.  CONCLUSIONS 


The  development  of  a  passive  floating  element  wall  shear 
stress  transducer  using  a  sensitive  capacitance-measuring 
technique  has  been  presented.  In  particular,  the  viability  of 
using  a  frequency  modulation  scheme  to  perform  the  capacitance 
measurement  has  been  demonstrated.  The  electronic  sensitivity  is 
such  that  the  transducer  is  able  to  measure  a  0.1  ym  displacement 
of  the  floating  element.  Expected  motions  of  the  floating  ele- 
ment are  less  than  5  ym  and  should  not  adversely  affect  the  char- 
acter of  the  boundary  layer  in  which  the  measurement  is  made. 
The  device  is  highly  linear,  enhancing  its  attractiveness  for  dy- 
namic wall  shear  stress  measurements. 


The  transducer  has  been  tested  and  found  to  predict  the 
steady  mean  wall  shear  stress  to  within  the  limits  of 
experimental  uncertainty  (66%);  two  of  the  three  data  points  ob- 
tained were  within  2  percent  of  the  actual  steady  mean  wall  shear 
stress.  The  transducer  must  be  calibrated  in  situ  in  order  to 
reduce  the  experimental  uncertainty;  workbench  calibration  is 
not  feasible  because  of  the  transducer's  exceptional  sensitivity. 
The  current  transducer  is  not  practical  for  dynamic  measurements 
of  the  mean  wall  shear  stress  since  it  does  not  possess  adequate 
bandwidth. 


The  gap  around  the  floating  element  has  been  verified  as 
being  the  greatest  single  vulnerability  of  the  floating  element 
transducer.  It  is  not  possible  to  seal  the  gap,  since  a  gap  seal 
in  the  current  design  immobilizes  the  floating  element  and  pre- 
vents transducer  response  to  wall  shear  stress.  On  the  other 
hand,  an  open  gap  affects  the  boundary  layer  and  introduces  er- 
rors of  its  own  to  the  wall  shear  stress  measurement. 
Fundamental  research  is  needed  in  this  area  in  order  to  fully 
understand  the  gap  effect. 


Finally,  a  course  of  action  has  been  recommended  to  improve 
the  capabilities  of  the  wall  shear  transducer  so  that  it  is  prac- 
tical for  dynamic  measures  of  the  shear  stress.  In  order  to  in- 
crease bandwidth,  reduce  transducer  sensitivity  to  inclination, 
and  minimize  the  adverse  effects  of  the  gap,  the  transducer 
should  be  microminiaturized.  Microelectronics  technology  is  be- 
coming increasingly  capable,  and  can  already  maintain  the  neces- 
sary sub-micron  tolerances. 
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APPENDIX  A 
CALCULATIONS 


A.l  Transducer  Model 

A. 2  Beam  Element  Equation  of  Motion 

A. 3  Transducer  Response 

A. 4  Transducer  Sensitivity  to  Inclination 
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A.l  Transducer  Model 


Using  the  model  of  Figure  A.l,  it  is  possible  to  develop  a 
general  theory  to  describe  the  output  of  the  wall  shear  transduc- 
er with  arbitrary  input.  Figure  A.l  highlights  the  important  fa- 
cets of  transducer  operation: 

-  the  wall  shear  stress  causes  a  deflection  of   the  beam 
element; 

-  deflection  of  the  beam  element  alters   its   capacitance 
with  respect  to  a  stationary  parallel  plate; 
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FIGURE  A. 1   Transducer  Model 
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the  change  in  capacitance  of  the  beam 
element-stationary  plate  changes  the  oscillation  fre- 
quency f  of  an  astable  multivibrator,  since  the  oscil- 
lation frequency  is  inversely  proportional  to  the  capa- 
citance; and 
-  the  frequency  difference  between  the  oscillator  signal 

and  a  reference  frequency  f     is  then  converted  to  a 

ref 

voltage  which  is  added  to  an  offset  voltage  V   ,  gener- 
ating the  transducer  output. 

The  capacitance  of  a  parallel  plate  capacitor  is  given  by 

C  =  K  ~  ^Hudson  (1944)1  (Al.l) 

where  K  is  a  constant  of  proportionality,  L  is  the  length  and  w 
the  width  of  each  plate,  and  s  is  the  distance  between  plates. 
If  one  plate  is  deflected  a  small  transverse  amount,  then  the 
distance  s  is  no  longer  constant  but  rather  a  function  of  x,  the 
coordinate  along  the  undeflected  length  of  the  movable  plate 
measured  from  the  wall.  Further,  the  deflection  may  also  be  a 
function  of  time.  Effectively,  s  not  a  constant  leads  to  an  in- 
finite number  of  parallel  capacitances 


dC  =  K   ,W  r  dx  (A1.2) 

s(x,t) 


in  the  limit  of  small  plate  deflections.   If  now  s(x,t)   is  de- 


79 


fined  as 


s(x,t)  =  dQ  +  y(x,t)  ,  (A1>3) 


where  d  is  the  undeflected  capacitor  plate  spacing  and  y(x,t)   is 
the  transverse  deflection  of  the  movable  plate,  then 

dC  =  K  1 — 7 ? TT  dx  =  K  — : ?'-\       •     f Al  U) 

d  +  y(x,t)        d  .   y(x,t)      ^Ai.**; 

d 


Nondimensionalizing  and  defining 


f  =   K  J  =  nU,t)  (A1.5) 


and  noting  that 


K-2—  =  C   ,  (A1.6) 

do 


where  C   is  the  capacitance  of  the  stationary-movable  plate  pair 
with  the  movable  plate  undeflected,  equation  Al.4  becomes 


dC  =  C   ;   ,d^  (A1#7) 


o  1  +  n(5,t)   * 


Equation  A1.7  is  valid  only  for  small  deflections  n(5,t)  ,  for 
which  three  dimensional  effects  are  insignificant.  The  total  ca- 
pacitance may  be  found  to  be 
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1 
C=   /dC  =Cof  l/na.t)  '  (A1'8) 


/' 


0 
"summing"  over  all  the  infinitesmal  capacitances  dC. 

The  oscillation  frequency  of  the  astable  multivibrator   is 
given  by 

f=Fc  (A1-9) 

where  C  is  the  capacitance  of  the  transducer  and  R  is  a  feedback, 
resistance  (constant).  Then,  for  a  deflection  of  the  movable 
plate  n(C,t),  the  oscillator  frequency  is 


,      .   1 

f=-L-  =  _L.  "l  (ALIO) 

R  C   RC 


J   i  +  n(C,t: 


or 


f  =  fn    C  HF  (ALU) 


O    f  dg 

J   i  +  n(?,t: 


Since 


V(t)  =V+G(f-f    )  (A1.12) 

o  ref  ' 


it  is  appropriate  to  select  f   _  =  f  to  obtain 

ref    o 


V 


(t)  =Vo  +  Gfo(f)"  >  (A1.13) 


or 


1 

V(t)  =  V  +  G  f      1  ,A1  ... 

o      o     /•    Jf  (A1.14) 

d? 


'/ 


1  +  n(C,t) 


using  equation  Al.ll. 


The  normalized  transducer  output  is  given  by 
V(t)  -  V, 


o 


G  f 
o 


t.    a 

y  i  +  n(?,t) 


-  1  •        (A1.15) 


0 


The  def lection n(?,t)is  determined  by  the  transducer's   excitation 
and  dynamics  considered  in  the  following  section. 


A. 2  Equation  of  Motion 


The  dynamics  of  the  movable  plate  of  the  transducer  deter- 
mines how  the  transducer  will  respond  to  applied  loads.  The  mov- 
able plate  is  cantilevered  at  its  base  and  has  a  floating  element 
mounted  on  its  free  end.  The  movable  plate  and  floating  element, 
together  called  the  beam  element,  may  be  modeled  as  in  Figure 
A. 2. 


F(t) 


M,  J 
y(x,t)  =  z(x,t)  -  z(0,t) 


•-  z(x,t) 


z(0,t)  =  z  (t) 


V/^.W/A 


FIGURE  A. 2   Beam  Element  Dynamic  Model 
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The  floating  element  is  modeled  as  a  mass  M  and  mass  moment 
of  inertia  J  lumped  at  the  free  end  of  the  cantilever  plate,  mo- 
deled as  a  uniform  beam  of  flexural  rigidity  EI  and  mass  per  unit 
length  y.  The  motion  of  the  beam  element  is  described  by  z(x,t) 
in  a  stationary  reference  frame.  In  order  to  develop  the  equa- 
tions of  motion  for  general  excitation,  a  concentrated  transverse 
force  F(t)  is  applied  to  the  floating  element  and  a  distributed 
transverse  force  f(x,t)  is  applied  along  the  beam. 


Of  interest  in  the  current  situation  is  the  motion  of  the 
beam  element  relative  to  its  base  rather  than  its  absolute  motion 
z(x,t),  since  the  stationary  capacitor  plate  is  actually  fixed 
with  respect  to  the  cantilever's  base.  In  effect,  it  is  desired 
to  obtain  the  equations  of  beam  element  motion  in  the 
non-inertial  reference  frame  translating  with  the  cantilever's 
base. 


If 


ZQ(t)  e  Z(0,t)  (A2.1) 


is  the  translation  of  the  base  of   the  beam  element,   then  the 
relative  motion  of  the  beam  with  respect  to  its  base  is 
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y(x,t)  =   z(x,t)  -  ZQ(t). 


(A2.2) 


Using  the  variational  technique  [Crandall,  et.  al.  (1968)]  and 
assuming  small  deflections  of  the  beam,  the  equation  of  motion  of 
the  beam  element  is  obtained  as 


d2Z 


"f£+«tE*+»dp-°  =  *<*.*> 


(A2.3) 


subject   to   the   boundary  conditions 

3, 


6y 


EI  tHj 

dX° 


=   0 


x  =   0 


92y 
d-xr 


(A2.4a) 
(A2.4b) 


=   0 


6y 


-   M 


•5 


d2: 


dt< 


33y 


+   EI 


&  +  ^ 


=    0 


(A2.4c) 


x  =   L 


J    3t23x 


-      EI 


3 


=    0 


(A2.4d) 


where    6  is   the  variational  operator.      For   the  case  where   the   beam 
is   cantilevered  at   its  base, 


<5y 


x=0 


(e) 


0 


(A2.5) 


x=0 


so  that  boundary  conditions  A2.4a  and  b  are   identically   satis- 
fied. 


03 


Clamped-Free  Beam. 

In  the  case  of  normal  transducer  operation,  the  beam  element 

3y 
is   free  at  x=L  so  that  the  variations  in  y  and  Tfr  at  x=L  do  not 

vanish  identically,  requiring  the  terms  in  brackets  in  equations 

A2.4c  and  d  to  be  equal  to  zero.   Thus, 


y  3 


s2    *      a  it  d2z 

^  ♦  EI  &     -      f  (,.t)  -  p  ^ 


(A2.6) 


with  boundary  conditions 

y(0,t)  =     0  (A2.7a) 

iZ(O.t)  ~=     0  (A2.7b) 

oX 

M -&»•'>  -  n  J&G-.O   =   F(t)   -  A  CA2.7C) 

3t"  3x°  dt^ 

_ETii2(L,t)      -   33y(L,t) 

EI^?       "  J  ~3F^3x"     '  <A2-7d> 


Note  that  acceleration  of   the   transducer  base   introduces 

D'Alembert   or   inertial   force  excitations   proportional   to 

d2z 

o  in  the  field  equation  (A2.6)  and  the  force  boundary  condi- 

'  It7    I 
tion   (A2.7c).    If   the   transducer  base   is  non-accelerating 

d2z 
i        -°  =  o  )>  the  transducer  is  kept  vertical  so  that  f(x,t)  =  0, 

and  the  floating  element  is  properly  aligned  so  that  there  are 


no 


error  forces  on  the  floating  element  (F(t)  =  x  (t)A  ) ,  then  the 

w    e 

field  equation  becomes 


oo 


y|^  +  EI  ^  -   0  (A2.8) 

subject  to  the  revised  boundary  condition  (A2.7c) 


Mi^a.t)   .  n>!jtt.t)   .  -  (A2.9) 

dt^  3x°  w     e 


with  the  other  three  boundary  conditions  unchanged.    Therefore, 

in  the  absence  of  the  error  terms  specified  above,  the  wall  shear 

stress  excites   the  homogeneous   response  of   the   beam  y(x,t) 

through  the  force  boundary  condition,  equation  A2.9.   Transverse 

acceleration  of   the   transducer  base  introduces  not   only  a 

d2z 

D'Alembert  force  at  the  floating  element  (x=L)  equal  to  -M  — t° 

j2  dt 

d^z 

but  also  a  distributed  D'Alembert  force  -y  -r-v     along  the   length 

dt 

of  the  beam.  Thus,  the  response  of  the  beam  element  to 
transverse  acceleration  of  its  base  is  dynamically  analogous  to 
its  response  to  the  wall  shear  stress  except  that  transverse  ac- 
celeration gives  rise  to  the  aforementioned  distributed 
D'Alembert  force. 


o/ 


Clamped-Plnned  Beam. 

If  the  beam  element  is  pinned  at  the  floating  element,  then 
fy  -  0  at  x=L  and  the  equations  describing  the  motion  of  the  beam 
element  are 

with  boundary  conditions 


y(o,t)  e  o 


(A2.7a) 


I**0'0   =   0  (A2.7b) 

dX 


y(L,t)  =  o  (A2'10) 


__  j£y(L,t)        d3y(L,t)  f 

~   EIlx^       "  JTFfc     '  (A2'7d) 


Because  F(t)  does  not  excite  the  transducer,  the  wall  shear 
stress  cannot  either.  Thus,  if  the  floating  element  exhibits  a 
pinned  end  condition,  the  beam  element  will  not  respond  to  the 
wall  shear  stress  although  it  will  respond  to  the  tranverse  ac- 
celeration of  the  transducer  base  (equation  A2.6  applies). 
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Clamped-Clamped  Beam. 

If  the  floating  element  exhibits  a  clamped   end  condition, 
then  equation  A2.6  is  still  the  appropriate  field  equation,  with 


Jz 

3x 
ment  will  not  respond  to  the  wall  shear  stress  but  will  respond 


y  and  — ^  identically  zero  at  x=0  and  x=L.   Again,  the  beam  ele- 


to  transverse  acceleration  of  the  transducer  base. 
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A. 3  Transducer  Response 


For  frequencies  much  smaller  than  the  fundamental  resonant 
frequency  of  the  beam  element,  the  deflection  shape  of  the  beam 
element  is  satisfactorily  described  by  its  static  deflection 
shape.  Assuming  that  the  transducer  base  is  non-accelerating, 
that  the  transducer  is  kept  vertical,  and  that  there  are  no  error 
forces  on  the  floating  element,  equation  A2.8  may  be  solved  with 
equations  A2.7a,b,  and  c  and  equation  A2.9  for  the  static  deflec- 
tion shape  (which  is  not  a  function  of  time).   Assuming 


F(t)   =  A  x  (t)  (A3.1) 

e  w 


leads  to 

A 
n(C,t)   =   2  k  ed  OK1   -   C3)  t  (t)  (A3. 2) 

s  o 

where  k  =  3El/L3is  the  effective  stiffness  of  the  beam  in   re- 

s 

sisting   an   applied   transverse   force   at  the  free  end. 
alternatively,  since 


A   _ 

£(t)  =  IT  Tw(t)  '    (A3-3> 

s 


is  the  transverse  deflection  of  the  free  end  of  the  beam, 
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nU.O  =  -±-  OK2     -  C3)  e(t) 
o 


(A3. 4) 


describes  the  deflection  shape  of  the  beam  element.   Since 


V(t)  -  V 


G  f 


1 


f        d6 

/  i  +  n(?,t: 


-   1    , 


(A1.15) 


0 


it  is  possible  to  evaluate  the  transducer  output  as  a  function  of 


\M  ■ 


A. 3.1  Linearized  Response 


If 


n(C,t) 


<<  1,  then 


/ 


d? 


1  +  n(?,t) 


-   -  l  +  /n(g,t)  d?  .   (A3. 5) 
£,t)  d£         0 


0 


Thus, 


V(t)  -  V 


G  f 


1 


t)  d£ 


(A3. 6) 


from  which,  with  equation  A3. 2, 


V(t)  -  V      A 
o      e 


G  f 


x  (t)  f 
:   d    J 


(3C2  -  53)dC  . 


2  k  d 

s   o   o 


(A3. 7) 


Upon  integration,  equation  A3. 7  leads  to 


yi 


V(t)  -  V       3  A 

«  —  t  (t)  (A3. 8) 

G  f         2  k  d   w 
o  so 


which  is  shown  in  Figure  2.5  as  a  function  of 


o       so 


the  normalized  end  deflection  of  the  beam, 


A. 3. 2  Exact  Response 

Using  equation  A3. 4,  the  integral  in  equation  Al.15   can  be 

d 
tegration  are  shown  in  Figure  2.5. 


E 

computed  numerically  for  different   j-  .   The  results  of  this  in~ 

o 
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A. 4  Sensitivity  to  Inclination 


When  the  transducer  is  inclined  from  the  vertical  by  an 
angle  9  ,  the  beam  element  deflects  by  its  own  weight.  For 
small  inclinations  from  the  vertical,  the  transverse  component 


beam  weight  per  unit 
length  w 


FIGURE  A. 3   Transducer  Inclination  Geometry 


F  of  a  load  F  is  given  by 

Ft  =  F  sin  i>   *   F  i>    ,  (A4.1) 

where  <P  is  the  local  angle  between  the  tangent  to   the  beam  and 
the  vertical  shown  in  Figure  A. 3.   Now, 


^  =  q    +  e  =  e  +  4^  (A4-2) 

o        o    dx 


dy 
where  9  =  -j   is  the  local  angle  between  the  tangent  to  the  beam 

and   the  undeflected  beam  position  and  9   is  the  angle  the  tran- 

r  o 

ducer  base  is  inclined  from  the  vertical. 


The  beam  element  consists  of  a  floating  element  of  mass  M 
and  a  cantilever  plate  of  mass  per  unit  length  u .  In  a  gravity 
field,  these  possess  weights  of  W  and  w,  respectively.  If  the 
transducer  is  inclined  from  the  vertical  with  no  shear  stress  ap- 
plied to  the  floating  element  and  without  accelerating  the  base 
of  the  transducer  transversely,  then,  from  equations  A2.6  and 
A2.7  and  by  considering  only  static  deflections  y(x)  of  the  beam 
element, 

EI  |^  =  f(x,t)  (A4.3) 

subject  to  the  boundary  conditions 


ya 


y(0)  =  0 


dy_(0)  = 


dx 


0 


-EI 


A(L>  =  F(t) 


dx^ 

d2y(L) 
dx2" 


=  0 


(A4.4a) 
(A4.4b) 
(A4.4c) 
(A4.4d) 


Both  f(x,t)  and  F(t)  are  developed  by  the  weight  of  the  canti- 
lever plate  and  the  floating  element  according  to  equations  A4.1 
and  A4.2: 


f(x,t)  =  w  (  e  +41  ) 

o    dx 


(A4.5) 


F(t)  -  W  (  6  +¥   ). 

o    dx 


(A4.6) 


Equations  A4.3  through  A4.6  are  then  solved  to  obtain 


n(0  =}Ci  +  C2  exp(6C)  + 


exp(-  -  O 


C:  sin 


/3 


3C  +  Cu  cos 


S3 


35 


-  ? 


e  l 

o 


(A4.7) 


wL 
where  9   is  in  radians  and  33  =  .   The  coefficients  C  were 

o  EI  i 

calculated  using  the  FORTRAN  program  GRAV3 ,  owing  to  the  complex- 
ity of  the  algebra  and  in  order  to  facilitate  the  ready  calcula- 
tion of  deflection  shapes  for  different  transducer  parameters. 
The  transducer  parameters  used  for  the  calculations  are: 

L  =  31  mm, 

dQ  =  375  ym  , 
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M  =  76  mg , 

M  ■  16  g/m,  and 


EI  -  2.84x10  Nm2 


for  which 


C  -  1.1568 
C  =  0.9807 

2 

C  -  2.1674 

3 

C  =  -2.1375 
4 

and  3  -  0.254314. 
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APPENDIX  B 
TRANSDUCER  NATURAL  FREQUENCIES 


B.l  Clamped-Free  Beam 

The  homogeneous  form  of  equations  A2.6  and  A2.7  may  be 
solved  in  order  to  determine  the  natural  vibration  frequencies  of 
the  clamped-free  beam  element: 


92v  .  „  94 


a 


-1+  EI  -fdf  =  0  (Bl.l) 

t         ot 


subject  to  the  boundary  conditions 

y(0,t)  =  0  (B1.2a) 

izCo.t)  3  (B1>2b) 

3x 


M 


i^a.t)  .     |^a,t)  (B1.2c) 


"§(l,l,»l^ll,-«  <»■"> 


Separating  variables  and  assuming  time  harmonic  behavior  leads  to 


y(x,t)  =  X(x)  T(t)  (B1.3) 


where 
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X(x)  =  Ci  sinh  Ax  +  C2  cosh  Ax  +  C3  sin  Ax  +  C4  cos  Ax   (B1.4a) 


and  T(t)  =  e"ia)t  (Bl.4b) 


with 


A^=^.  (B1.5) 

EI 


Imposing  the  boundary  conditions  (Bl.2)  obtains 


cosh  AL  +  cos  AL  -  y(sinh  AL  -  sin  AL)      M   .   _  _ 

sinh  AL  -  sin  AL  -  y(cosh  AL   -  cos  AL)      y  (Bl.ba) 


where 


w2J 


sinh  AL  +  sin  AL  -   ..    (cosh  AL  -  cos  AL) 


Y  = 


AH .       (B1.6b) 


U)2J 

cosh  AL  +  cos  AL  -  -r-rrr  (sinh  AL  +  sin  AL) 

AE1 


The  roots  of  equations  Bl.6  determine  the  eigenvalues  A. 
The  first  three  natural  frequencies  have  been  calculated  and  are 
tabulated  in  Table  V.  The  effect  of  the  floating  element's  mass 
moment  of  inertia  was  approximated  by  modeling  the  element  as  a 
slab  4  mm  by  4  mm  and  2  mm  thick,  and  was  found  to  be  negligible 
for  the  first  three  vibration  modes. 
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B.2.   Clamped-Pinned  Beam 

The  determination  of  the  resonant  frequencies  of  the 
clamped-pinned  beam  follows  the  procedure  for  the  clamped-free 
beam  except  that  equation  Bl.2c  is  replaced  by 

y(L,t)  =  0   .  (B2.1) 

The  characteristic  equation  then  becomes 

sinh  XL  -  sin  XL   -  y(cosh  XL  -  cos  XL)  =  0       (B2.2) 

where  y  is  as  given  in  equation  Bl.6b.  The  natural  frequencies 
of  the  first  three  clamped-pinned  modes  is  tabulated  in  Table  V, 
both  including  and  neglecting  the  effect  of  the  rotational  iner- 
tia of  the  floating  element.  Blevins'  (1979)  tabulation  of  the 
natural  frequencies  for  the  clamped-  pinned  beam  are  the  same  as 
the  natural  frequencies  predicted  from  equation  B2.2  above. 


B.3  Clamped-Clamped  Beam 

For  the  clamped-clamped  case,  the  mass  moment  of  inertia  no 
longer  enters  the  problem  since 


iy_(L,t)  5  0  (B3.1) 

3x 


is  the  new  boundary  condition  that  replaces  equation  Bl.2d.  The 
natural  frequencies  of  the  first  three  clamped-clamped  modes  are 
taken  from  Blevins  (1979)  and  shown  in  Table  V. 
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TABLE  V 


Summary  of  Beam  Element  Natural  Frequencies 


Clamped-Free  Beam      f 

f2 

f3 

J  =  0           61.7 

414.2 

1202 

J^O           61.7 

411.2 

1180 

A%                0 
Clamped-Pinned  Beam    f. 

0.7% 
f2 

1.9% 
f3 

J  =  0          340.6 

1104 

2303 

J  4   0          339.4 

1091 

2246 

A%             0.4% 
Clamped-Clamped  Beam   f 

1.1% 
f2 

2.5% 
f3 

494.2 

1362 

2671 

Notes: 


1. 
2. 

3. 

4. 


All  frequencies  are  in  Hz. 

-4     2 
J  4   0  indicates  J  assumed  to  1.267x10   mg  m  , 

modeling  the  floating  element  as  a  rectangular  slab 

4  mm  x  4  mm  x  2  mm. 

The  floating  element  mass  M  and  (where  appropriate) 
mass  moment  of  inertia  J  are  modeled  as  being  lumped 
at  the  undamped  end  of  the  beam  element. 

Clamped- clamped  data  and  verification  of  clamped-pinned 
calculations  are  from  Blevins  (1979). 


APPENDIX  C 
UNCERTAINTY  ANALYSIS 


From  equation  2.1, 


8  k   d 

Tw  S  5— ^  (V  "  V  '  (C  1) 

W    3  A  f  G         °  lC,JJ 

e  o 


It  is  shown  in  Cook  and  Rabinowicz  (1963)  that  if 

q  =  F(a±)  (C.2) 


where  F  is  a  function  of  n  independent  variables  a.  ,   then   the 

2 
most  probable  mean  square  error  e   in  the  quantity  q  is 


•  t  (■£-)'  <  ■ 

i-1  V  9ai7     X 


ea  =   Z^  I- "J   e;.  .  (C3) 


assuming  the  a.  are  statistically  independent.  Each  term 
ec  denotes  the  error  or  uncertainty  in  the  quantity  C .  From 
equation  C.3, 


e 
x  a 


where  e.=  —  and  the  uncertainty  in  V(t)  has  been  neglected  since 


it  is  an  order  of  magnitude  smaller  than  e.T 

vo 


3wt3 


Since  k  =  3EI/L0  =  TT7~  and  A  =  wz  ,  equation  C.3  obtains 


4L 


LUZ 


£2     -     £2     +     £*     +      (3e     )2     +      (3e     )2 

E  w  t  L  (C.5a) 

-    (.05)2    +    (.05)2    +    9(.01)2    +    9(.05)2    -.0284      (c.5b) 


and 

ek   =    2z{    =    2(.05)2    =    0.0050  (c6) 

so  that,  from  equation  C.4, 


zl    a    0.0284    +    (.05)z    +    0.005    +    (.01)2    +    (.05)2    + 
(.10)2[ ±£Jl }2    «    0.429 


and 


/ 2^^_\ 

\(1   Pa)  (.4   V/Pa)/ 


e      =    0.66 

T  • 


(C.7) 


(C.8) 


n 
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